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ABSTRACT: New models of the Fock space sector corresponding to some 
fixed number of electrons are introduced. These models originate from the rep- 
resentability theory and their practical implementation may lead to essential re- 
duction of dimensions of intermediate CI spaces. A certain zero-order theory that 
gives wavefunctions approximately equivalent to ones obtained by accounting all 
excitations from the Hartree-Fock reference state up to the g-th order is proposed. 
Simple numerical examples are given to illustrate our approach. 
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1. Introduction 

In present work new realizations of p-electron sector of the Fock space 
are proposed and their simplest properties are studied. Notations and ter- 
minilogy of this article are new for Quantum Chemistry and require therefore 
few preliminary comments. The basic idea of our approach is widespread in 
modern mathematics: to assemble information about an object as a whole 
using available local information. This idea was properly formalized and 
systematically used by Jean Leray who introduced his famous sheaves over 
topological spaces and employed them in complex analysis (see, e.g., [0]). 
Since this branch of mathematics is not used in Quantum Chemistry , at 
least at present, we found it possible to borrow existing mathematical terms 
"sheaf" and "germ" filling them with quite different meaning. Namely, we in- 
troduce a family of g-electron functions (germs) and call such family "sheaf" 
if the germs constituting this sheaf are subject to certain "gluing" conditions. 
Then it is shown how sheaves thus defined can be used to restore unique p- 
electron wavefunction. The whole theory was motivated by the theorem from 
the pure representability theory that was proved in . To make presentation 
of our approach more transparent and logical, in the second section necessary 
basic definitions and results from the representability theory are formulated. 
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Our approach is essentially finite-dimensional and set-theoretical manipula- 
tions and technique of enumerative combinatorics are extensively used. For 
these reasons we found it possible to change notations that seems to be con- 
sidered as traditional in literature on the represent ability problem. Namely, 
capital letters are normally used to denote subsets of the spin-orbital index 
set N and/or the orbital index set M. Number of elements (indices) in index 
sets is usually denoted by lower case letters, e.g. |A^| = n, \M\ = m, etc. 
With such notations it seems consistent to use lower case letters for the num- 
ber of electrons and for the current density matrix order that are just the 
numbers of elements in the relevant spin-orbital index sets. Throughout this 
article n is the number of molecular spin orbitals (MSO), m is the number of 
molecular orbitals (MO), p stands for the number of electrons, and q is the 
density matrix order. 

In the third section the formal theory of (p, g)-sheaves is presented. 

In the fourth section construction of correlated basis functions in config- 
uration interaction (CI) theory is discussed. 

In the fifth section it is demonstrated how the abstract mathematical 
theory of {p, g)-sheaves may be used for replacing large-scale CI calculations 
by series of CI calculations in relatively small dimensions. 
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2. Basic Definitions. 



Let J^ri,i be one-electron Fock space spanned by an orthonormal set of n 
molecular spin-orbitals . Electronic Fock space is defined as 



— ^^n,p (1) 

p=0 



where 



p 

^ n,p = /\^n,l (2) 
^n,0 = C (3) 

and C is the field of complex numbers. 

"Determinant" basis vectors of the Fock space will be labelled by subsets 
of the spin-orbital index set N: for any R <Z N the corresponding basis de- 
terminant will be denoted by \R). Note that basis determinants are labelled 
by subsets and all sign conventions connected with their actual representa- 
tion as the Grassman product of ordered spin-orbitals are included in the 
definition of the creation- annihilation operators: 

al\R) = {l-6,^R){-ir\RUt) (4a) 
a,\R) = 6,,R{-ir\R\t) (46) 

where 

^ ri, ifiei? 

may be considered as a possible generalization of the Kronecker S symbol, 
and 

e= |{l,2,...,i-l}ni?| (6) 

is the sign counter. 

Dirac's ket-bra reahzation of the operator space over the electronic Fock 
space {Endc{J-n-i ^n) ~ ® ^n) equipped with the trace inner product 

{z\t) = Tr{zH) (7) 

will be used. 
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The contraction operator over JF„ (g) JF* is defined as 

n 

c = ^ a| (8) 

i=l 

Its restriction to the p-electron sector of the operator space is identical 
(up to the nonessential combinatorial pref actor) to the commonly used con- 
traction operator ||, |, 0|. 

Let tp G J-'n,p ® ^n,p be some p-electron operator. It can be expanded via 
basis "determinant generators" as 

(p) 

tp = Y.tRs\R){S\ (9) 

R,S 

where the upper summation index (p) indicates that the sum is taken 
over all p-element subsets of the spin-orbital index set A^. It can be shown 
0, § that 

(p) (fe) 

c\tp) = k\Y.tRs E (10) 

R,S KcRns 

The set-theoretical operation Ak was introduced in ^ and its main prop- 
erties were investigated in [0] but due to its importance for the uniform treat- 
ment of the sign prefactors arising in the expressions of the type of Eq.(lO) 
we repeat here necessary definitions. 

For two arbitrary index sets R and 5* their symmetric difference is RAS = 
{RU S)\{Rn S) . This binary operation on the set of all subsets of some given 
set is commutative and associative with the empty set as its unit. For any 
subset K = {ki, k2, ■ ■ ■ , ki} C N we put 

Ax = {l,2,...,fci}A{l,2,...,A;2}A...A{l,2,...,fcJ. (11) 

For example, if K = {2, 4, 5} in = {1, 2, 3, 4, 5, . . .} then 

Ak = {1,2}A{1,2,3,4}A{1,2,3,4,5} = {3, 4}A{1, 2, 3, 4, 5} = {1,2,5} 

The electronic Hamiltonian associated with some chosen spin-orbital basis 
set in one-electron Fock space is of the form 

H = E(^l^li)4«i + ^ E {ij\ — \kl)ala]aiak (12) 

i,j ^ i,j,k,l ^12 
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The energy functional on the Fock space is defined as 



E{t) = Tr{Ht) (13) 

and is a hnear mapping from J-'n ® J-'n to the field of complex numbers. 

Using specific form of the electronic Hamiltonian, it is possible to contract 
the energy domain and redefine the energy in terms of g-electron operators 
(g > 2). Introducing reduced Hamiltonian that acts on the g-electron sector 
of the Fock space 

Hp^g = - — I I](«|/i|j>laj + ^ H {ij\ — \kl)ala]aiak (14) 

^ i,j ^ i,j,k,l ^"12 

we can recast the energy expression (13) in the form 





E{tp) = ^Tr 



(15) 



where tp is some j9-electron operator. On the right hand side of this equation 
g-electron image (2 < g < min{p,n—p)) of tp with respect to the contraction 
appears. 

Of actual interest in Quantum Chemistry are the so-called pure p-electron 
states that are associated with p-electron wavef unctions. The corresponding 
pure representability problem was formulated in the very first articles of 
Coleman i i i 0: 

(PR) To find the necessary and sufEcient conditions that should be im- 
posed on q-electron operator tg to guarantee that there exists at least one 
p-electron wave function \1/ such that |ycP~'^|\l/)(\l/| = tg. 

From the very beginning it was understood that the structure of the set 
of all pure representable g-electron operators is very complicated and that 
its constructive description that could be of any practical value hardly exists 
at all. The well-known theorem of the Convex Sets Theory states that a 
real-valued linear function defined on some compact convex set necessarily 
reaches its minimal value at one of the extreme points of this function domain. 
Therefore, if the set of all pure representable g-electron operators is replaced 
by its convex hull, the set of certainly more simple structure will be obtained 
and one can guarantee that the direct search of energy minimal value on this 
set will lead to the optimal pure state. This was also employed by Coleman 
who formulated his famous ensemble representability problem: 
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(ER) To find the necessary and sufEcient conditions that should be im- 
posed on q-electron operator tq to guarantee that there exists ensemble of 
p-electron states such that ^c^~'^J2^i\'^i){^i\ = tq- 

Here I] = 1 and Aj > 0. 

i 

It is very easy to find three very simple conditions that should necessarily 
be satisfied by any representable g-electron operator: 

(i) Nonnegativeness; 

(ii) Hermiteancy; 

(iii) Fixed trace value (usually Tr(tq) = 1). 

The set of all g-electron operators satisfying these three conditions will 
be denoted by £n,q- Operators from £n,q are called density operators. 

In his pioneer works Coleman [0, ^, ^ found the solution of the ensemble 
represent ability problem for the case q = 1 and formulated two very strong 
necessary conditions for the case q = 2. 

In our previous work we generalized part of Coleman's results to treat 
the case of density operators of arbitrary order and found the following de- 
scription for the set Vn,p,q giving an exterior approximation of the set of all 
ensemble representable density operators of order q: 

V„,p,<j = A'^{n,p, q)£n,q n £n,q (16) 

where the automorphism A{n, p, q) is determined by its matrix representation 
as 

\q-s)\ q ) K'CN\{IUJ) [\KnK'\) 

(17) 

with respect to specially selected operator basis 

e[^ = (-l)l(^uJ)nA^||jU^^^jU^| (^g) 

Here / n J = 0, s = |/| = \ J\, and ajj = |(/ U J) fl A7v\(/uJ)|- In the same 
work very simple connection between the operator A{n,p,q) and its inverse 
was found: 

A~^(n,p, q) = A{n,n — p, q) (19) 
Let us apply the operator A{n,p,q)^cP~'^ to some pure state |\I')(^E'| where 

(p) 

vl> = ^ Cn\R) (20) 
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Rather comphcated manipulations (see Appendix D of 0) lead to 

1 1 (P+9) 

A(n,p,g)|c^-'^|vl/)(vl/| = -py E d'liZ) (21) 



ZcN 



where 
and 



dqiZ) = \i;z)m (22) 

|^z)= E(-l)l(^\")"^-IC^z\5|5) (23) 
scz 



Let us introduce the set 

Bn,p,q = {{Z, S) C N X N : \Z\ = p + qL\S\ = qLS C Z} (24) 
and the equivalence relation on this set 

{Z, S) ~ {Z', S') ^ Z\S = Z'\S' (25) 

The set of all equivalence classes Bn-p^q contains exactly elements and in 

each equivalence class (Z, S) there are ("~^) elements. 

Eqs.(21)-(23) together with the definition (24) make it possible to formu- 
late the following necessary and sufficient conditions of pure represent ability. 
Theorem . q-electron operator tq is representable by pure p-electron state 
if and only if its image dq witli respect to A{n,p,q) satisfies the following 
conditions: 

(i) Tliere exists expansion 

I (P+9) 

d, = 7;-^i: dq{Z) (26) 



ZcN 



such that for each Z G N, dq{Z) is either the null operator or corresponds 
to (unnormalized) pure q-electron state \ipz); 
(a) The mapping 

{Z, S) ^ (_l)l(AS)nAzl ^^1^^^ (27) 



is constant on the equivalence classes {Z, S). 
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Unfortunately, condition (i) of the above Theorem is depressingly non- 
constructive. Only for the case p + q = n conditions (i) and (ii) are trivial 
and the following corollary of the main theorem may be formulated. 

Corollary 1. For p + q = n, q-electron operator tg is representable by 
pure p-electron state if and only if its image dq with respect to A{n,p,q) 
corresponds to pure q-electron state. 

The next corollary shows how to construct p-electron wavefunction from 
q-electron ones. 

Corollary 2 If the conditions of Theorem are fuUfilled then (up to nor- 
malization ) the required pure state may be presented as 

\^) = ^{-l)\^^\'^'^'''\{S\^z)\Z\S) (28) 

where the sum runs over equivalence classes of the set Bn,p,q modulo the 
equivalence relation (25). 

The above Theorem was first formulated and proved in our previous arti- 
cle 0. Its Corollary 1 has analogue in the ensemble representability theory 
where the constructive description of the set of all representable density op- 
erators of order q for the case n = p + q was well-known for years |^ . 

At first glance the above Theorem is hardly of any practical value since 
in general case we do not have even the slightest idea how to expand den- 
sity operator dq. In the next section it will be shown that it is a wrong 
impression and that this Theorem opens a completely new direction both in 
mathematical and computational chemistry. 
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3. Representation of Many Electron Wave Function As a 
Sheaf of q-Electron Germs. 

For any (p + g)-element subset Z of the spin-orbital index set N we put 

(g) 

^„,,(Z) = 0C|S) (29) 

SCLZ 

Definition. Family {V'zjzcAf of g-electron functions, il)z G J^n,q{.Z), is called 
a sheaf of q-electron germs of p-electron wavefunction, or just a (p, g)-sheaf, 
if the mapping 

(Z, S) ^ (-l)l(^\^)"^-l(5|^z), \Z\=p + q, \S\ = q,SGZ (30) 

is constant on the equivalence classes of the set -Bn,^,^ modulo the equivalence 
relation (25). 

The set of all (p, g)-sheaves will be denoted as Sn,p,q- 

The Theorem formulated in the preceding section and its Corollary 2 
lead to the conclusion that there exists one-to-one correspondence between 
p-electron wavefunctions and {p, g)-sheaves for each fixed q. Indeed, ar- 
bitrary wavefunction \Ef represented by its Full CI (FCI) expansion (20) is 
corresponded with {p, g)-sheaf 

ZCN (31) 

where g-electron wavefunctions ipz are given by Eq.(23). On the other hand, 
due to the conditions of gluing (30) each (p, g)-sheaf uniquely determines p- 
electron wavefunction in accordance to Eq.(28). The mapping Sn,p,q performs 
disassembling of p-electron wavefunction into the family of g-electron ones 
whereas its inverse assembles p-electron function from its g-electron germs. 

Disassembling mapping (31) may be used to transfer all relevant struc- 
tures from the Fock space J-'n^p to the set Sn,p,q- The most important linear 
structure is transfered as 

Sn,p,,(Ec.^O = |(Ec.*Oz| =T.(^d^z}zcN (32) 
i |- i J ZcN i 

^By an abuse of terminology arbitrary vector of p-electron sector of the Fock space is 
called wave function even if it is not normalized. 
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whereas the inner product (Euclidean structure) is defined as 

I ip+q) 

{{Mzcn\{(I)z}zcn) = 7— Y {'^z\(t>z) (33) 

U j 

To avoid cumbersome notations, we agree to omit, when possible, subscript 
Z C in writing elements from Sn,p,q- 

With the aid of the disassembling mapping (31) it is possible for each 
Gxed q to identify the Fock space jF^^p with the space of all (p, g)-sheaves 
^n,p,q considering the later just as an isomorphic model of the FCI space. 

The mapping 

'^n,p,g{Z) : {ipz'}z'cN ^ V'z (34) 

is obviously linear and performs a projection of the vector space Sn,p,q on 
the vector space J^n,q{Z). Among its sections (right inverses) there exists in 
a certain sense canonical one that is defined as follows. To each g-electron 
wavefunction 

(9) (p) 
V'z = E Cs\S) = E (-l)l^"^^l^ii|^\i?) (35) 

SCZ RcZ 

where R — Z\S and 

C'«=(-l)l«nA.|^^ (36) 
it is possible to put into correspondence the unique ]9-electron wavefunction 

(p) 

*z - E Cr\R) (37) 

RcZ 

with its subsequent disassembling. The resulting section of the projection 
(34) will be denoted by jn,p,q{Z) and is a hnear injective mapping satisfying 
the standard relation 

T^n,p,qiZ) O jn,p,q{Z) = idyr^ ^(^z) (38) 

where idj^^^i^z) is the identity operator over J^n,q{Z)- As follows from this 
equation, the image of a certain g-electron wavefunction ipz with respect to 
jn,p,q{Z), that will be denoted as {ipzz'}z'cN , posesses the property : 

^z = ^zz (39) 
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or, in the other words, the initial g-electron wavefunction ijjz is necessarily a 
germ of the (p, g)-sheaf {ipzz'}z'cN- This property motivates the following 
definition. 

Definition, (p, q)-sheaf {ipzjzcN is called simple if it is generated by 
one of its germ, that is, if there exists {p + q)- element subset Z' C N such 
that 

{MzcN^jn,pAZ'){M (40) 

Proposition 1. {p, q)-sheaf generated by any its nonzero germ corresponds 

to single determinant wavefunction. 

Proof. If sheaf {ipz}zcN is generated by any its germ then the wave- 
function (37) does not depend on concrete subset Z. If we suppose that it 
involves different determinants \R) and \R') in its expansion, then it is pos- 
sible to find {p + r/)-element subsets Z, Z' G N that separate it! and R', that 
is R C Z, R' C Z',hut R Z' and K <f. Z. In this case, however, ^z can 
not be equal to ^^'B 

Definition, (p, q)-sheaves generated by any its nonzero germ are called 
determinant sheaves. 

It is clear that determinant sheaves constitute an orthonormal basis of 
the vector space Sn,p,q- Selection of this basis will lead to the approach 
completely equivalent to the FCI one in the basis of p-electron determinants. 
All possible advantages and disadvantages of our model for FCI space may 
appear in this case either on the matrix element evaluation level or on the 
level of algorithms for solving partial eigenvalue problem. Concrete formulas 
for matrix elements ready for practical implementation of the conventional 
CI method are collected in Appendix A. Here more general approach that is 
not based on the determinant sheaves is discussed. 

g-electron density operator dg of the form of Eq.(26) is associated with 
each (p, g)-sheaf {ipz}zcN- Recalling that dg originates from some g-electron 
density operator tg = A~^{n,p,q)dq and that the operator A{n,p,q) is Her- 
mitean in trace inner product, we can recast general energy expression (15) 
in the form 

£;({^z}) = Tr [{A-\n,p,q)Hp^g) dg = 
1 (p+«) 

J—. E {i^z\A-\n,p,q)H,^g\i;z) (41) 

[ g ) ZCN 

Note that in this formula the initial definition of the reduced Hamiltonian is 
slightly modified and the current operator Hp^g differs from that in Eq.(14) 
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(p) 

by the combinatorial pref actor jf / _ 

In exactly the same manner general expression for the matrix elements 
of the Hamiltonian between two (p, gf)-sheaves may be written 

1 (p+g) 

H{{M,{<Pz}) = 7— Y E {i^z\A-\n,p,q)H,^,\<f>z) (42) 

[g ) ZCN 

This completely general formula expresses arbitrary matrix element of p- 
electron Hamiltonian via matrix elements of the transformed q-electron re- 
duced Hamiltonian. It is necessary to emphasize that Eq.(42) does not couple 
g-electron spaces J-'n,q{Z) with different Z. 

Of course, for their practical use it is necessary to rewrite Eqs.(41)-(42) 
in terms of molecular orbitals taking into account simplest spin symmetry 
restrictions. Relevant formulas are collected in Appendix A. 

Let us consider some arbitrary (p, g)-sheaf {ipz}zcN- Each its nonzero 
germ ^/^z^. generates simple {p, g)-sheaf {ijJZjzjzcN = 3n,p,q{Zj){ipZj)- The set 
of all such simple sheaves will be denoted as Jn,p,q{{'^z}z<zN) or Jn,p,q{^) 
where ^ is the wave function corresponding to the initial (p, (3')-sheaf. The 
subspace of the vector space Sn,p,q generated by the set Jn,p,q{^) will be called 
the g-subspace associated with the wave function ^> (or the corresponding 
(p, q')-sheaf) and will be denoted as IV„,p,g(^): 

Wn,p,qm=Y.^{i^Z,z}zcN (43) 
3 

It is pertinent to note that the sum on the right-hand side of Eq. (43) is not 
direct. 

The importance of the definition (43) is due to the following fact. 
Proposition 2. {i/jzjzcN e Wn,p,q{^) for any (p, q)-sheaf {ipz}zcN = 

Sn,p,g(^)- 

Proof. Eq.(28) may be recast as 
1 (p+i) (p) 

l^)= A^E E(-i)"''^^'(^W^)l^) 

\ q ) ZCN RCZ 

Prom Eqs.(35)-(37) it follows that 

Sn,p,q I E(-l)l^"^-l(Ai?l^z)l^) I = Jn,p,q{Z){i^z) 
\RCZ J 
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But jn,p,q{Z){iJz) is, by definition, a simple (p, g)-sheaf generated by ipz- As 
a result, s„p,j(^) is a linear combination of its simple sheaves and, conse- 
quently, belongs to VFn,p,g(^) ■ 

Thus, to each ]9-electron wavefunction it is possible to put into correspon- 
dence the family of subspaces 

* - (W^n,p,l(*), Wn,,,2{n W^n,p,,_(*)) (44) 

and the associated integral vector 

^ ^ind{^) = {tndi{^),tnd2{^), . . . ,tndg^^^{^)) (45) 
where Sn,p,q is the disassembling mapping (31), Qmax — min{p,n —p), and 

indq{^!) = dimcWn,p,q{'if) = rankcJn,p,q{^) (46) 

Definition. indg{^) is called the q-index of the wavefunction ^ (or the 
corresponding (p, q)-shcaf). 

At present we do not know if indices (46) can be considered as integral 
characteristics of the associated wave function. It is certainly true for the 
so-called formal Fock space. 

Definition. Vector space of all formal linear combinations of all subsets 
of the index set N with complex coefficients is called formal Fock space. 

Formal Fock space is in a certain sense universal model of the electronic 

Fock space with orbital specific erased. In this space vector (45) is the vector 

of integral chracteristics of formal p-electron wavefunction ^. 

On the Fock space additional structure connected with the action of the 

p 

group of unitary transformations /\u induced by one-electron ones appears. 
The same wavefunction has, in general, different CI coefficients in different 
MSO bases and, consequently, is corresponded with different sheaves. If the 
q- subspaces associated with these sheaves can be of different dimensions is 
the question to be investigated. 

When q is fixed, then the index (46) may be called just the CI index 
of the wave function under consideration. The meaning of this index is the 
following: for arbitrary wave function the length of its expansion in the 
subspacc Wn,p,q{'^) is not greater than indq{^) (for fixed MSO basis set). 

Let us consider the group Gn,p,q{^) of automorphisms of the vector space 
W^n,p,g(^) leaving invariant the finite set of simple sheaves J„^p^g(^). Since 
Jn,p,q{'^) generates H^n,p,q(^) this group is isomorphic to the subgroup of the 
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symmetric group permuting simple sheaves from J„ p^^(^), and, consequently, 
is finite. If for a given wave function ^ there exist non-trivial groups Gn^p^qi"^) 
is the problem for future study. At present we can state only that for each 
q one-dimensional subspace of Wn^p^qi"^) spanned by Sn,p,q{'^) carries trivial 
representation of Gn,p,q{^) (see Eq.(28)). 

Definition, (p, q)-sheaf {i'z} = Sn,p,q{^) is called stable if it coincides 
with the lowest eigensheaf of the CI matrix in the subspace Wn^p^q{'^). 

It is easy to see that any sheaf can be transformed to the stable one in 
no more than one step. 

Definition. Wavefunction ^ is called stable if its (p, q)-sheaves Sn,p,q{^) 
are stable for each 2 < q < min{p, n — p). 

All determinant functions and FCI lowest eigenvector are necessarily sta- 
ble. In general case with each j9-electron wavefunction from JF„ p we can 
associate the unique family of its stable images ^q from Sn,p,q {q > 2) and 
the following inequalities always keep true:(^|if|\l') > {'^q\H\'^q). 
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4. Construction of Correlated g-Electron Basis Functions. 



It is well-known that the proper selection of many electron basis func- 
tions becomes the crucial stage of the CI approach when the FCI limit is 
unreachable. The commonly used methods of basis functions construction 
employ very robust physical idea of successive generation of k-fold excitations 
from some reference state (/c = 1,2, . . .). In the simplest case the Hartree- 
Fock (HF) single determinant is used as a reference state. We start our 
analysis with the conventional CI approach in the basis of the determinant 
(p, g)-sheaves. 

Selecting in the subspace J^n,q{Z) the function 

and applying the sheaf generation scheme described by Eqs.(35)-(37) lead to 
the determinant (p, g)-sheaf 



iZ',S')^{i-l)^^'^'^\\Z'\R)} (47) 



where R = Z\S, and it is immediately clear that for any pair {Z', S') G {Z, S) 
the same determinant sheaf will be obtained. This result can be reformulated 
in the following manner: determinant {p, g)-sheaves are in one-to-one corre- 
spondence with the equivalence classes of the set Bn,p,q modulo the equiv- 
alence relation (25). Selecting HF determinant as the reference state, it is 
easy to construct basis determinant sheaves by applying Eq.(47) succesively 
to the "excited pairs" {Z, S) ^ Z\S = R\I U J. This at first glance trivial 
scheme can be generalized in a very interesting manner. 

Let us try to select q'-electron wavefunctions ipz £ J'n,q{Z) with the fol- 
lowing property: the corresponding simple sheaf should account as much of 
correlation effects as possible. To this end we should turn to Corollary 1 of 
our Theorem from the second section. From this corollary it follows that for 
fixed Z the electronic energy minimum is reached on the sheaf generated by 
the lowest eigenvector of the operator 

A-\p + q,p,q){VzH,^,Vz) (48) 

where Vz '■ ^n,q ^n,q{Z) is the standard projection. 

Moving along some equivalence class {Z,S), one can select the lowest 
eigenvector ip(z,s) of the operator (48) obeying if necessary some additional 
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conditions of the type {i>(z,s)\'5) 7^ 0. Then all simple (p, gf)-sheaves |'0(z,5)z'| 
are to be constructed and the lowest eigenvector of the CI matrix in the vec- 
tor space generated by these sheaves should be found. The resulting sheaf 
should be made stable by applying the routine described in preceding sec- 
tion. With such an approach to each equivalence class {Z, S) some vector 
accounting part of the correlation effects is put into correspondence. These 
(in general, non-orthogonal) vectors are in one-to-one correspondence with 
p-electron determinants but the problem of their linear independence is still 
to be studied. 

Applying this scheme to the equivalence class of the HF p-electron de- 
terminant leads to a certain zero-order theory corresponding approximately 
to the conventional CI method accounting excitations up to the q-th order. 
The main difference between the standard CI approach and our scheme is in 
replacing large CI matrix partial diagonalization by a series of partial diago- 
nalizations of relatively small orders. Indeed, for class representatives matri- 
ces of one of the following orders arise: (^'^q-k^)-> where k — 0,1, . . . ,q. 
After finding optimal vector for each class representative, it is necessary to 
solve partial diagonalization problem of no more that order to get 
equivalence class energy and wavefunction. And the last diagonalization to 
stabilize the equivalence class wavefunction is performed in the space of the 
minimal possible for a given wavefunction dimension equal to its CI index. 
Numerical tests are discussed in detail in the next section. Here simple com- 
parison of relevant dimensions is performed. For the fixed value Ms of the 
total spin projection the number of a electrons Pa and f3 electrons Pf^ is fixed 
and the total number of all excitations up to the q-th order from the HF 
reference state in space of m molecular orbitals is equal to 




For the case of 80 electrons correlated in the one-electron space of 200 MOs, 
q — 2, and the total spin projection Ms — there are Z(i^2) = 60816000 
excited configurations. At the same time for each representative of the HF 
determinant equivalence class partial eigenvalue problem of order either 1681 
or 861 should be handled. To get the wavefunction of the HF equivalence 
class it is necessary to find the lowest eigenvalue of the matrix of order not 
larger than 51040. 

To improve the approximation described one can try to account relevant 
excitations from the zero-order wavefunction. We believe, however, that it 



16 



is possible to develop very efficient new algorithms based on our model of 
the CI space. Indeed, for some current (p, g)-sheaf {ipz}zcN approximate 
Newton-Raphson (NR) shift vector can be calculated with its subsequent 
disassembling it to give NR sheaf {4>z}zcn- Then in the vector space gener- 
ated by simple sheaves {ipzz'}z'cN and {(f)zz'}z'cN the CI problem should 
be solved to to give new stable sheaf, etc. With such an approach both p- 
and g-electron functions are actually optimized and it can be expected that 
the final FCI solution will be obtained by subsequent passing through CI 
subspaces of moderate sizes. 
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5. Simple Numerical Examples. 



To illustrate our approach, we performed very simple calculations of the 
HF equivalence class wavefunctions for two small molecules in ST0-6G Gaus- 
sian basis set. There are two reasons why such primitive examples were 
chosen. Firstly, the whole approach is just at its very beginning and sepa- 
rate probably very complicated and time consuming work should be done to 
develop efficient computer schemes and create the corresponding computer 
code. Secondly, for small molecules it is easy to write down explicitely all 
necessary equivalence classes, their energies, etc., to make understanding of 
previous sections of this article easier to the reader who has only vague im- 
pression about equivalence classes and is not experienced in set-theoretical 
manipulations. In all calculations Handy's split determinant representation 
was used 0. In this representation spin-orbital index set is replaced by a pair 
of orbital index sets for a and /3 spins. Binary codes for index sets identifi- 
cation were also employed: e.g., index set {1,3,5} as a subset of the orbital 
index set M = {1, 2, 3, 4, 5, 6, 7} is represented as a binary code 1010100. All 
RHF calculations to generate MO basis set were performed with the aid of 
the GAMESS program 0. 

From the definition (24) of the set -Bn.p.g and the equivalence relation 
(25) it is easy to see that the same (p + g) -element subset Z C N occurs in 
exactly (^^''^) different equivalence classes. At the same time for fixed single 
determinant reference state \R) its equivalence class representatives can be 
uniquely labelled by subsets Z D R, or, in orbital representation, by pairs 
{Za, Zfs) D (Ra, Rfs) such that \Za\ + iZ/sl = p + q. 

For the LiH ground state {RliH = 1.595 A) in the basis of 6 MOs the 
HF code is \HF) = (110000, 110000) and its equivalence class for g = 2 
contains exactly (^^2~^) ~ elements. These elements {Za, Zp) are of three 
different types: type (1,1) of pairs \Za\ = \Zp\ = 3, type (2,0) of pairs 
\Zo,\ = 4, iZpl = 2, and type (0,2) of pairs \Za\ = 2, |Z^| = 4. For each 
representative of the HF equivalence class we selected the lowest eigenstate 
of the operator (48) as a germ generating simple (4,2)-sheaf. Of 28 simple 
sheaves only 9 proved to be linearly independent (most probably due to high 
point symmetry) and the equivalence class energy was calculated to give E = 
—7.972047 a.u. The corresponding wavefunction is not stable and after its 
disassembling 35 linearly independent simple sheaves were generated to give 
finally stable wavefunction of the CI index 35 with energy E = —7.972323 
a.u. Conventional approach leads to CI space of 92 excited determinants and 
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the energy value E = —7.972323 a.u. for the lowest CI vector in this space. 
For the water molecule ground state (Ron = 0.957 A. HOH = 104.3°) 

HF) = (1111100, 1111100). 
6 elements of types (1,1), (2,0), and 



in the basis of 7 MOs the HF reference code is 

/14-10 



Its equivalence class contains 
(0,2). The corresponding codes and energies are: 



Pair (Z^, Zp) 



1111110 
1111110 
1111110 

1111101 
1111101 
1111110 



Energy 
of HF class 
representative 
(a.u.) 

-75.693408 
-75.687226 
-75.687226 



Energy 
of HF class before 
stabilization 
(a.u.) 



Energy 
of HF class after 
stabilization 
(a.u.) 



-75.716895 



-75.728024 



-75.695390 



-75.680388 



-75.680388 



1111101 
1111101 
1111111 
1111100 
1111100 

1111111 

The CI wavefunction involing HF determinant and all single and dou- 
ble excitations from this determinant as a reference state (on the whole 141 
determinants) corresponds to the energy E — —75.728063 a.u. In our calcu- 
lations we selected the lowest 2-electron wavefunctions for each representative 
of the HF class by diagonalizing operator (48): for (1,1) type elements the 
order of the eigenvalue problem was 36 (four matrices), for (2,0) and (0,2) 
type elements we had two matrices of order 21. All six simple (10,2)-sheaves 
proved to be linearly independent and after solving 6 x 6 CI problem we 
obtained wavefunction (sheaf) of the CI index 45. Final stabihzation of this 
wavefunction was performed in the CI space of dimension 45. 

In Tables 1-2 integral characteristics of ground state wave functions of 
LiH and H2O molecules obtained by CI method accounting all single and 
double excitations from the HF reference state sd) and by FCI method 
are presented. In these tables term 'actual expansion length' stands for the 
numbers of expansion coefficients greater than 10~^ by absolute value. 



19 



The examples discussed are certainly not very impressing and their only 
advantage is in their simplicity. But even on these primitive examples it was 
possible to learn what computational difficulties may occur in the course of 
our approach implementation to treat larger systems. The ffist difficulty is 
connected with binary codes (index subsets) generation. Standard generation 
in the lexical order is known to be slow |[Tl|. Much faster Gray algorithms 
| rT| should certainly be employed in future implementations. Technique to 
determine code number if the binary code is given, and restore code by its 
number in the list of codes (without scanning the list) is also required. The 
second difficulty is the selection of linearly independent simple sheaves and 
their orthogonalization. The problem here may be connected with large 
amount of required simple sheaves and their numerical orthogonalization. 
An attempt should be made to develop some analytic technique for selection 
of linearly independent simple sheaves. And at last it is necessary to have 
fast routine for the Hamiltonian matrix elements evaluation in the basis of 
sheaves obtained by orthogonalization of simple sheaves. 
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6. Conclusion 

Our approach to treatment of correlation effects based on new models of 
the FCl spaces opens new perspective directions both in mathematical and 
computational chemistry. Staying on the zero-level theory (HF equivalence 
class) it is possible to improve the quality of the wavefunctions by increasing 
the order of density matrices moving from {p, 2)-sheaves to sheaves of higher 
order. On the other hand, it is possible to perform calculations in the Fock 
space model Sn,p,q with fixed q and try to reach the FCI limit either by turn- 
ing to the 1st, 2nd and higher order theories or use from the very beginning 
purely algorithmic approach close in its idea to the Davidson diagonaliza- 



tion scheme |jT2|]. There are many things to be done in the frameworks of our 
approach both in pure theory and in practical implementation of new compu- 
tational schemes. Sheaves properties such as their point and spin symmetry, 
definition of excited sheaves, etc., are still to be studied. Of special interest 
are {p, l)-sheaves and their possible connection with the density functional 
theory approach. Vector and integral characteristics of formal many electron 
wavefunctions obtained in this article are very intriguing as well. 
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Appendix A. 



Operator A^^{n,p,q) involved in general expression (42) for the CI ma- 
trix elements is defined by its matrix representation with respect to the basis 
of g-electron determinant generators (18) and due to the huge dimension of 
space where it acts it may be very difficult to perform its direct transforma- 
tion to some other, probably more convenient, basis. Therefore, at present 
stage of theory, all required matrix elements should finally be expressed via 
matrix elements of Hamiltonian between determinant {p, (3')-sheaves. For two 
determinant sheaves (45) corresponding to determinant labels R and R' in 
spin-orbital representation we have 

HnR' = j^ (-l)l^^^-l+l«'"^-l(Z\i?|A-^(n,p,g)if,^,|Z\i?') 

\ q ) ZdRUR' 

(Al) 

In the majority of cases all molecular calculations are performed in one- 
electron bases of molecular orbitals. Turning to MO basis in Eq.(A.l) leads to 
more cumbersome expressions but simplifies essentially concrete calculations. 

Using Handy's split representation of spin-orbital index sets [|], it is pos- 
sible to replace R by (i?^, Rjs) and Z by {Za, Zp) where R^ and Z„ are taken 
from the same orbital index set M, \M\ = m for cr = a, (3. With fixed projec- 
tion Ms of the total spin the numbers Pa of « and p/j of [3 electrons are fixed 
for each basis determinant and the following simple equalities keep true: 

\Z^\R^\ + \Zp\Rp\ = q {A.2) 

\Z^\ + \Zp\=p + q (A3) 

Let us start with the analysis of the sign prefactors in Eq.(A.l). To this 
end it is convenient to introduce the following notations: 

= K\{R^ n K) = {n,r,, . . . ,rfcj (AA) 

K = R.\iR. n K) = {si, S2, . . . , SkJ {A.5) 

where \K„\ = \K'^\ = k„. 

Technique for handling sign prefactors developed in our previous article 
makes it easy to get the following comparison modulo 2: 

|i?n A^l + |i?'n A;.| = 
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= |(i?«Ai?:,)nA^J + |(i?^Ai?^)nA^J + |(i?,Ai?^| xmod(|Z^|,2) (A6) 
Let us put 

where X^- is a "free" part of the summation index Z^- in Eq.(A.l) and 

\X^\ + \Xf,\=q-{ka + k0) (A8) 

Now the sign prefactor power can be separated into two parts: Non-active 
part that is not involved in summation 



E 



\RaAR'^\ + l 



+ J2\{R^AK)nAn^nK 



(A9) 



and active part 

a = ^ I {RaAR'^) n AxJ + \RaAR'„\ x mod(|i?;3 UR'pUXp\,2) (A. 10) 

Brackets in Eq.(A.9) stand for the entier function extracting integral part of 
its argument. 

The comparision 

|i?^A7?;| = 2A;^ = O(mod 2) 
leads to the following expressions that hardly admit further simplification 

(All) 



a 

a = ^|(i?,Ai?;)nAxJ 



(A12) 



Now it is possible to rewrite general expression (A.l) for the CI matrix 
elements in terms of orbital index sets: 



:-ii 



\ q J '«''/3 Xo,CM\{RcUR'J 

(ia+'/3=q-('=a+fc/3)) 



^ i-iy{K^UX^,KpUXp\A-\2m,p,q)Hp^,\KuX^,K'f,UXfs) 

Xf,cM\iRi3UR'^) 

{A.13) 
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Note that for fixed spin projection value both the reduced Hamiltonian Hp^q 
and its image A^^{n,p, q)Hp^g are block-diagonal with exactly q + 1 blocks. 

Eq.(A.13) can be simplified for each concrete pair {ka, kp). In the most 
important case q = 2 we have six such pairs:(0,0),(l,0),(0,l), (1,1), (2,0), and 
(0,2). 

I. (A;a,M = (0>0) 

This case corresponds to — R'^ and 



_ 1 

{Ra,R/3),{Ra,R/3) ^2m-p 



^ {ij,$\A-\2m,p,2)H,^2\im+ 



i,jeM\Ra 



+ {^,tj\A-\2m,p,2)H,^2\^,tj) + 



(A14) 



+ E E (i,i|^-'(2m,p,2)i/,_2|i,i) 

ieM\Ra jeM\Rf} 

II. {k^,k^) = (1,0) 

In this case = {ri},K'^ = {si},ri si, £ = 1 + \{ri, Si} n Ar^^r'J, 



i}, 0|A-'(2m,p, 2)//p_2|{si, i}, 0) + 

ieM\{RaVJR'^) 



and 



+ ^ (ri,i|^-i(2m,p,2)i/p^2|si/ 

i&M\Ri3 



0, if ri < i and si < i or ri > i and si> i 



(^.15) 
(A16) 



1, if ri < i and Si > i or ri > i and s\ < i 
III. (fc„,A;^) = (0,l) 

In this case Kf) = {ri},K'^ = {si},ri si, £ = 1 + |{ri,Si} n A^^^nR^I, 



2 J 



^(Ra ,Rg ) , (Ra ,R'g ) 

Y (-1)"(0, {n, i}\A-\2m,p, 2)H,^2% {si, 0) + 

ieM\{RpLiR') 
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+ J2 {t,n\A'\'2m,p,2)Hp^2\i,Si) , (A17) 

iGM\Ra 

and a is given by Eq.(A.16) 

IV. iK,k^) = ii,i) 

We have = {ri},K'^ = {si},ri =^ si, Kp = {rs},/^^ = {s2},r2 ^ S2, 
^ = \{rusi} n AR^nR'J + |{'^2,S2} H Ar^h/?;,!, a = 0, and 

H(R^,R^),(R'^,^^) = j^^^{ri,r2\A \2m,p,2)Hp^2\si, S2) (A18) 

V. {K,kfs) = (2,0) 

In this case Ka = {ri,r2},K'^ = {si,S2}, e = |{ri, rs, Si, S2} n A^^^nR^I, 
a = 0, and 

HiR^,R,),iR:^,R,) = p^{nr2,^\A-\2m,p,2)Hp^2\siS2,^) (A19) 
V. ik^,kfs) = (0,2) 

We have Kp = {ri,r2},K^ = {si,S2}, £ = |{ri, r2, Si, S2} n A^^ni?;, I, 
a = 0, and 

= p^{(l>,nr2\A-\2m,p,2)H,^2\(l>,siS2) (A20) 

Practicaly all large scale CI calculations are based on the so-called direct 
methods and use Davidson-type algorithms for solving partial eigenvalue 
problem [l^. In these methods at each Davidson iteration the result of 
applying Hamiltonian to the current CI vectors is calculated and explicit 
Hamiltonian construction is actually not required. In the frameworks of our 
approach it is not difficult to write down general formula for Hamiltonian 
(12) action on some trial (p, g)-sheaf {ijjz} (in MSG basis). We have 

{R\m = j;^^j:{-ir''^^KZ\R\A-\n,p,q)H,^,iJz (A21) 
zdR 



n—p 
<? 



where \R) is some basis determinant in CI wavefunction expansion and V'z are 
g-germs of obtained by its disassembling. This formula, after its rewriting 
in MO basis set, can certainly be used in existing direct CI algorithms but 
we believe that it is possible to develop very efficient new methods based on 
our model of the CI space. One of such algorithms is briefly outlined in the 
fourth section. 
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Table 1: LiH ground state wave functions: Integral characteristics 

















q 


Total 


CI 


Actual 


Total 


CI 


Actual 


value 


number 


index 


expansion 


number 


index 


expansion 




of germs 




length 


of germs 




length 


1 


208 


35 


25 


360 


69 


57 


2 


482 


35 


25 


720 


69 


52 


3 


584 


35 


24 


752 


69 


47 


4 


436 


35 


10 


495 


69 


15 



Table 2: H2O ground state wave functions: Integral characteristics 







'^sd 






^FCI 




q 


Total 


CI 


Actual 


Total 


CI 


Actual 


value 


number 


index 


expansion 


number 


index 


expansion 




of germs 




length 


of germs 




length 


1 


128 


49 


32 


232 


120 


67 


2 


69 


45 


14 


91 


61 


21 


3 


14 


12 


7 


14 


12 


7 


4 


1 


1 


1 


1 


1 


1 
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